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We report the results of Chandra and XMM-Newton observations of an ultra-luminous supersoft X-ray source 
in M101. M101 ULX-1 underwent 2 outbursts in 2004 during which the peak bolometric luminosities reached 
10 41 ergss -1 . The outburst spectra were very soft and can generally be fitted with a blackbody model with 
temperatures of 50-160 eV. In two of the observations, absorption edges at 0.33 keV, 0.56 keV, 0.66 keV, and 
0.88 keV were found. A cool accretion disk was also found in the 2004 December outburst. During the low 
luminosity state, a power-law tail was seen up to 7 keV. It is clear the source changed from a low/hard state 
to a high/soft state. In addition, it showed at least 5 outbursts between 1996 and 2004. This is the first ultra- 
luminous X-ray source for which recurrent outbursts with state transitions similar to Galactic X-ray binaries 
have been observed. From the Hubble Space Telescope data, we found an optical counterpart to the source. 
During the 2004 outbursts, we also performed radio and ground-based optical observations. All the results 
strongly suggest that the accreting object is a > 28 OOM 0 black hole. 


1. Introduction 

Recent high angular resolution X-ray observations 
reveal that there are a large number of ultra-luminous 
X-ray sources (ULXs) in many nearby galaxies. ULXs 
are luminous (Lx > 10 39 ergss -1 ) non-nuclear X-ray 
point sources with apparent X-ray luminosities above 
the Eddington limit for a ~ 10 Mq black hole (BH). 
While some ULXs have been associated with super¬ 
novae, many are thought to be accreting objects with 
X-ray flux variability observed on timescales of hours 
to years. A natural possibility is that the compact ob¬ 
ject is an intermediate-mass black hole (IMBH) with 
mass of ~ 1O 2-4 M 0 i- The origin of such objects 
remains uncertain. Some ULXs may have a stellar- 
mass black hole with beamed emission 0! While 
current observations are inconclusive about the nature 
of ULXs, recent X-ray observations show that some 
ULXs have a cool accretion disk (kT ~ 0.1 keV), sug¬ 
gesting the presence of IMBHs 3J. 

While the majority of ULXs have X-ray emission 
from 0.1 to 10 keV, a few ULXs have vervsoft spectra 
with no X-ray emission above 1 keV 0013; similar to 
supersoft X-ray sources (SSSs) in the Milky Way. The 
high luminosities of ultra-luminous SSSs are inconsis¬ 
tent with typical nuclear burning white dwarf models 
for Galactic SSSs. These ultra-luminous SSSs could, 
however, very well be IMBHs. Their luminosities and 
temperatures are consistent with what is predicted for 
accreting BHs with masses between roughly 100 and 
1000 Mq. Alternatively, outflows from stellar-mass 
BHs could also achieve such a high luminosity and 
low temperature @]. 


In this paper, we report a series of X- 
ray/optical/radio observations of the ultra-luminous 
SSS in M101 (CXOU J140332.3+542103; M101 ULX- 
1 hereafter) during the low states and outbursts. 


2. X-ray Observations 

M101 ULX-1 is one of the most luminous ULXs. It 
was discovered with ROSAT and was confirmed as a 
SSS with a blackbody temperature of about 100 eV, 
with Chandra 0 Eft Ell ■ During 2000 March, Chan¬ 
dra detected it at Lx ~ 4 x 10 39 ergs s -1 , and then in 
2000 October, its luminosity was around 10 39 ergs s -1 . 
In 2004, Chandra conducted a monitoring program 
for M101. Figure 1 (Left) shows the long-term X-ray 
lightcurve of M101 ULX-1 from 2004 January to 2005 
January. M101 ULX-1 was near the detection limit 
during January, March, and May; the X-ray spec¬ 
tra were harder with a power-law shape (see Figure 1 
Right), and the X-ray luminosity was about 2 x 10 37 
ergs s -1 , a factor of about 10 2 fainter than that dur¬ 
ing the outbursts in 2000 E2- The source was found 
to be in outburst during the July 5 observation, with 
an X-ray luminosity of about 7 x 10 39 ergs s -1 . Data 
taken on July 6, 7, and 8 show that the source was 
in outburst with a peak bolometric luminosity (for as¬ 
sumed isotropic emission) of about 10 41 ergs s -1 j!3 |. 
In general, the X-ray spectra are best described with 
an absorbed blackbody model with temperatures of 
~ 50—100 eV (see Figure 1). In addition, we found ab¬ 
sorption edges at 0.33, 0.57, 0.66, and 0.88 keV in two 
of the high state spectra. These features may signal 
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Figure 1: Left: The light curve of M101 ULX-1 in 2004 as observed by Chandra (solid circles) and XMM-Newton (open 
circles). It is very clear that the source stays at ~ 10 37 ergs s -1 during the low luminosity state. When it went into an 
outburst in 2004 July and December, the peak 0.3-7 keV luminosity was near 10 41 ergss -1 , a factor of ~ 1000 
comparing to the low state. We also show the optical (INT and CFHT) and radio (VLA) coverage in the figure. It is 
worth noting that the CFHT observation was taken simultaneously with XMM-Newton. Right: Unfolded spectra of 
M101 ULX-1. The total spectrum, blackbody component, and power-law component are shown in green, red, and blue, 
respectively. During the 2004 July outburst, the spectra can be fitted with blackbody/disk blackbody model with 
temperatures of 50-100 eV. In two of the observations, absorption edges at 0.33, 0.57 (July 8), 0.66, and 0.88 keV (July 
6) were found. On July 7, the peak bolometric luminosity reached 10 41 ergss -1 (which is one of the most luminous 
ULXs). The XMM-Newton spectrum taken on July 23 can be fitted with a blackbody (kT = 53 eV) and a power-law 
(a = 0.72) model. The combined Jan-May spectrum is similar to the XMM-Newton spectrum, but with much lower 
luminosities (~ 10 37 ergss -1 ). 


the presence of highly ionized gas in the vicinity of the 
accretor (e.g., warm absorber). A DDT XMM-Newton 
observation was made on July 23 and the luminosity 
was about 6 x 10 38 ergs s -1 . A harder X-ray spectrum 
with a power-law tail (kT = 53 eV and a = 0.72) 
was seen up to 7 keV. More recently, Chandra obser¬ 
vations made in September and November indicated 
that the source returned to the low state with a power- 
law spectrum and a luminosity of ~ 2 x 10 37 ergs s -1 . 
In 2004 December, the source was in outburst again 
with very soft spectra (kT = 50 — 160 eV) and a peak 
luminosity of about 10 41 ergs s -1 . During the rise of 
the outburst, the spectra were supersoft with black¬ 
body temperatures of 40-70 eV. The source showed 
a cool accretion disk near the peak of the outburst 
on 2005 January 1. The X-ray spectrum can be fit¬ 
ted with a disk blackbody model with kT^BB = 0.16 
keV which is significantly harder than previous super- 
soft X-ray spectra. The peak luminosity of the 2004 
December outburst is about 9 x 10 40 ergss -1 which 
is very similar to the 2004 July outburst. A DDT 
XMM-Newton observation made on 2005 January 8 
indicated that the source was in the decline stage 
(Lx = 6 x 10 39 ergss -1 ) and the X-ray s pect rum re¬ 
turned to a supersoft state (kT = 56 eV) llil . 


3. Optical Observations 

M101 ULX-1 is located near star forming regions 
in a spiral arm. Figure 2 shows a far UV (1350- 
1750 Angstroms) image of M101 taken with GALEX. 
While GALEX cannot detect M101 ULX-1 because of 
the sensitivity and confusion, the ultra-luminous SSS 
is about 10" away from a bright UV source. Based 
on the position provided by Chandra, we searched for 
optical counterpart of M101 ULX-1 in archival Hub¬ 
ble Space Telescope (HST) data. The region of M101 
ULX-1 was observed with HST using Wide Field Plan¬ 
etary Camera 2 (WFPC2) and Advanced Camera for 
Surveys (ACS) in 1994, 1995, and 2002. After reduc¬ 
ing the data with standard procedures and correcting 
the astrometry of the HST and Chandra images, we 
found a blue object (V = 23.8, B — V = —0.2) within 
the 0.6" Chandra error circle (Figure 3). At the dis¬ 
tance of M101 (d=6.7 Mpc), the absolute magnitude 
corresponds to My = —5.3. The magnitudes and col¬ 
ors of the blue star are consistent with an OB star. 
During the 2004 July outburst, we performed a se¬ 
ries of target-of-opportunity observations at ground- 
based optical telescopes (see Figure 1 and 3). These 
observations included WIYN 3.5m, the 4.2m William 
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Figure 2: GALEX far UV (1350-1750 Angstroms) image 
of M101. The blue circle is the location of M101 ULX-1. 


Herschel Telescope (WHT), the 2.5m Isaac Newton 
Telescope (INT), and the 3.6m Canada-France-Hawaii 
Telescope (CFHT). Unfortunately, WIYN and WHT 
observations were affected by bad weather and bright 
moon. M101 ULX-1 was barely detected in the INT 
images but was clear seen in the CFHT images even 
under very bad seeing (Figure 3). It is worth noting 
that the CFHT observations were performed simulta¬ 
neously with XMM-Newton. These are the first simul¬ 
taneous X-ray/optical observations for an ULX. Dur¬ 
ing the 2004 December outburst, we also performed 
an optical follow-up observation with the WIYN 3.5m 
on 2005 January 17 and M101 ULX-1 was clearly de¬ 
tected. 


4. Radio Observations 

After the discovery of the outburst of M101 ULX-1 
in 2004 July with Chandra, we observed M101 ULX- 
1 under target-of-opportunity time with Very Large 
Array (VLA) D-array at 4.86 (14 arcsec resolution) 
and 8.46 GHz (8.4 arcsec resolution; Fig. 4) for about 
1.5 hours on 2004 July 21. There is no radio coun¬ 
terpart with a 3cr detection limit of 0.075 mJy/beam 
at 8.46 GHz. We also searched for the radio emission 
from archival VLA data. In the 2002 May 7 A-array 
1.4 GHz (1.5 arcsec resolution) data, the source was 
not detected in the 6-hour observation but there is a 
2(7 hint of a source about 1" away from M101 ULX- 
1; the 3cr detection limit of the observation is 0.075 
mJy/beam. 


5. Discussion 

M101 ULX-1 is a very unique ULX. It has a very 
soft blackbody spectrum during outburst while it has 
a hard power-law tail in the low state. It also has 


recurrent outbursts which provide an excellent oppor¬ 
tunity to explore a totally new variability timescale for 
an ULX. Only a handful of transient ULXs have been 
detected. This may be due to the lack of repeated 
observations of nearby galaxies. M101 ULX-1 is the 
only ULX that has been observed regularly with X- 
ray telescopes. Hence, rapid multiwavelength follow¬ 
up observations become possible and meaningful. The 
X-ray/optical/radio observaitons of M101 ULX-1 re¬ 
ported here is the first and only rapid follow-up ob¬ 
servations for an ULX. 

Although over 100 ULXs have been discovered Llj, 
only a handful of sources have Lx> 10 41 ergs s _1 , 
assuming isotropic radiation [lj, [l]i [lQ, [l7|. A lu¬ 
minosity this high is difficult to achieve in an X-ray 
binary unless the accretor has a mass greater than 
roughly 10 M 0 . We do have evidence that M101 ULX- 
1 is likely to be an X-ray binary, since its luminosity 
has been observed to change by a factor of > 10 on a 
time scale of hours 0,0 While it is possible that 
our luminosity estimates are higher than the true lu¬ 
minosity, the effects that lead to overestimates, such 
as beaming, or various anisotropies, tend to change es¬ 
timated luminosities downward by a factor of roughly 
10. Its unusually high bolometric luminosity (10 41 ergs 
s~ 4 ), coupled with its short-time-scale time variabil¬ 
ity, therefore make M101 ULX-1 a good candidate for 
an accreting IMBH. For instance, anisotropic X-ray 
emission can result super-Eddington luminosity for 
a stellar-mass BH 0. However, in order to achieve 
such a high luminosity for a stellar-mass BH, extreme 
beaming is required and the disk is expected to be 
much hotter. Similarly, the temperature of radiation 
pressure-dominated accretion disk model proposed is 
too high [l3|. The pure blackbody spectrum also 
makes relativistically beaming unlikely .3j. 

Its soft spectrum in the high state is another im¬ 
portant piece of evidence. An accretion disk around 
a very massive BH is expected to produce supersoft 
X-ray emission |u|. If we used the 90% lower lim¬ 
its of the inner disk temperature derived from the 
disk blackbody fits, the BH mass is estimated to be 
> 28OOM 0 |2£|. This is consistent with prior work 
on IMBH models for ULXs 0. It also complements 
other work on evidence for cool disks in ULXs that 
has also been considered as evidence for IMBH mod¬ 
els 01 although a power-law component usually con¬ 
tributes a significant fraction of X-ray emission. In¬ 
deed, M101 ULX-1 showed a cool accretion disk spec¬ 
trum (kTpBB = 0.16 keV) during the December out¬ 
burst. Furthermore, the high state luminosity is ap¬ 
proximately 16% of the Eddington luminosity for a 
6000 Mq BH; we therefore expect the inner disk to 
be optically thick, which is consistent with the IMBH 
interpretation. 

The state changes of M101 ULX-1 span a luminos¬ 
ity range larger than a factor of 1000. M101 ULX- 
1 was first detected by ROSAT in 1996 November 
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Figure 3: Upper Left: HST /ACS B band image (900s) of the field near the SSS. The image has been corrected for 
astrometry using 2MASS catalog. The 0.6” Chandra error circle is indicated. The B magnitude is about 23.6. Lower 
Left: HST/WFPC2 U band (F336W) image (1200s) of the same field. The U magnitude of the source is about 22.3. 
Upper Right: CFHT g' band image (600s) of the same field. The image was taken simultaneously with XMM-Newton 
on 2004 July 23. Lower Right: CFHT v! band image (600s) of the same field. 


[2lj . with an extrapolated unabsorbed 0.3-7 keV lu¬ 
minosity of 8 x 10 38 ergss -1 (blackbody model with 
Nr = 10 21 cm -2 , kT = 75 eV). Higher absorp¬ 
tion (3 x 10 21 cm -2 ) would imply a luminosity up 
of 4 x 10 39 ergss -1 . The source was not detected in 
other ROSAT observations |2lj. In a 108 ks ROSAT 
HRI observation taken in 1996 May/June, we can set 
a 3 ct detection limit at 6.8 x 10 37 ergss -1 (power-law 
model with Njj = 1.2 x 10 2 ° cm -2 and a = 2). Setting 
the Nh at 1.5 x 10 21 cm -2 , the luminosity limit be¬ 
comes 1.5 x 10 38 ergss -1 . In 2000 March and October, 
Chandra detected the SSS in the very soft state with 
luminosities of ~ 10 39 ergss -1 0, El I n 2002 June, 
XMM-Newton observed M101 and did not detect the 
SSS [22]. We re-analysed the image and found that the 
SSS may be barely detected with possible contamina¬ 
tion from nearby sources. We derived the 0.3-7 keV 
luminosity of the SSS as ~ 10 3 ' ergss -1 , by assuming 
a power-law model. The source was also in the low 
state between 2004 January and May with luminosity 
of ~ 10 37 ergss -1 . It is clear that the SSS has had at 
least 5 major outbursts (Lx > 10 39 ergss -1 ). On the 
other hand, the very low luminosities (~ 10 37 ergs s -1 ) 


during 2002 and 2004 indicate that the source was in 
the low state. The source varies by as much as a fac¬ 
tor of 1000 between the low state and the high state. 
This amplitude is even greater than many Galactic 
BHs. Remarkably, the source also shows spectral 
changes. It is very clear that there is a power-law 
component (a = 1.4) in the composite low state spec¬ 
trum, while the high state spectra are supersoft. The 
XMM-Newton observation taken during the decay of 
the July outburst revealed that the blackbody com¬ 
ponent was still strong and there was a very hard 
power-law tail, similar to the low state. In addition, 
the source showed spectral change during the 2004 De¬ 
cember outburst. During the peak of the outburst, the 
X-ray spectrum chan ged from a supersoft spectrum to 
a quasisoft spectrum 1211 . It is clear that the spectrum 
was significantly harder than typical outburst super- 
soft spectra. This is the first and only ULX for which 
recurrent outbursts with state transitions similar to 
Galactic X-ray binaries have been observed. 

With a peak bolometric luminosity ss 10 41 ergss -1 , 
it is likely that IMBH is the central engine of the sys¬ 
tem. Interestingly, the blackbody component is al- 
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Figure 4: VLA observations of M101 ULX-1 with 4.86 GHz (14 arcsec resolution; Left) and 1.4 GHz (1.5 arcsec 
resolution; Right). The cross is the Chandra position of M101 ULX-1. 


ways seen while the power-law component becomes 
stronger when the source is at lower luminosities. It 
may indicate that the power-law component is due to 
Comptonization of soft photons [24}. The photon in¬ 
dex is, however, harder than that of Galactic BHs in 
the low state (a ~ 1.7). None of the Galactic BHs has 
similar spectrum [25}. The closest example is Galac¬ 
tic microquasar V4641 Sgr for which the photon index 
was measured to be between 0.6 and 1.3 1251261. 

Another important feature of our X-ray spectral fits 
is the presence of absorption edges, which we found at 
0.33, 0.57, 0.66, and 0.88 keV in two of the high state 
spectra. These are consistent with C V, N VI, N VII, 
and O VIII edges. We note, however, that the 0.33 
keV edge may be due to calibration of the ACIS-S 
near the carbon edge around 0.28 keV. These features 
may signal the presence of highly ionized gas in the 
vicinity of the accretor (e.g., warm absorber), and may 
be consistentwith an outflow from the source. In fact, 
outflow models have been suggested for ULXs in which 
the soft component of some ULXs is due to outflow 
of a stellar-mass BH Q. Similar argument to explain 
the super-Eddington luminosity of M101 ULX-1 [11} . 
However, the new data which reveal the changes of 
temperature and bolometric luminosity, the extremely 
high luminosity, and the state transition of the source 
will be difficult to explain by such a model. It re¬ 
mains a puzzle that we saw different edges in the two 
observations, possibly related to the geometry of the 
system. Nevertheless, outflows from IMBHs may be 
expected. 

It is also interesting to note that similar absorption 
edges are expected and have been observed in white 
dwarf (WDs) systems such as the recurrent nova U 


Sco |3, CAL 87 and RX J0925.7-4758 |H|. SSS ra¬ 
diation with luminosities in the range from 10 3 ' ergs 
s^ 1 to roughly the Eddington limit for a 1.4M@ WD 
are expected for nuclear burning WDs 0. Indeed, 
such models have been proposed to explain SSSs in 
the Galaxy and Magellanic Clouds. Note, however, 
that neutron star and BH models have also been pro¬ 
posed. The high state luminosity of this source is al¬ 
most two orders of magnitude larger than that of any 
other known SSS, and it rules out steady nuclear burn¬ 
ing on hot WD models as an explanation for M101 
ULX-1. A similar argument can be made for neutron 
star models. A nova explosion could explain some fea¬ 
tures of the data. It is not favored, however, because 
the photospheric radius and lack of hard radiation at 
the peak of the outburst have no obvious explanation. 

We expect that at least some of the ULXs involving 
IMBH accretors are transients [ 10 - The required 
condition is that the donor must be a massive star (> 
5 M 0 ) in regions of young populations. M101 ULX- 
1 satisfies these conditions. The colors of the optical 
counterpart are consistent with an OB star. More 
recent optical observations suggest that the optical 
counterpart is a B supergiant with a mass of 9 — 12M 0 
and the optical spectrum indicates that M101 ULX-1 
is consistent with a high-mass X-ray binary l3ll . The 
source is also very close to star forming regions in a 
spiral arm as indicated in Figure 2. 

It may be impossible to conclusively establish, using 
current technology, that any X-ray source in an exter¬ 
nal galaxy is an IMBH. A system like M101 ULX-1 is 
therefore particularly valuable, because it provides 5 
different pieces of evidence that together make a con¬ 
sistent argument in favor of an IMBH interpretation. 
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Its high-state luminosity, its short-time-scale variabil¬ 
ity, its soft high-state spectrum, its pronounced spec¬ 
tral changes, and the high-mass companion all suggest 
that M101 ULX-1 is a strong IMBH candidate. 
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